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Abstract: 

A new mechanism for hard inclusive production, which leads to a violation of the fac- 
torization theorem, is suggested. The mechanism is illustrated by a detailed discussion of 
Higgs meson production in high energy deutron-deutron scattering. Numerical estimates 
^ • f° r the effect are given for high energy hadron (nucleus) scattering. 
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In this letter we consider a new mechanism for hard inclusive production, which vio- 
lates the factorization theorem [ij. We illustrate this mechanism for the case of inclusive 
heavy Higgs production in deutron-deutron high energy scattering. The suggested mech- 
anism is present in any central hard inclusive production, such as the high inclusive 
production of mini jets and jets. 

The usual description of inclusive production of a Higgs meson in a deutron-deutron 
interaction, is given by a simple Mueller diagram M shown in Fig. la. Our calculation is 
based on the impulse approximation in which the production mechanism is initiated by a 
primary single nucleon-nucleon collision. Assuming the factorization theorem, the cross 
section of interest is given by 

a(Higgs) = J dx 1 dx 2 F£(x 1 , M 2 )Fg(x 2 , M 2 ) a(hard) , (1) 

where F^(x, M 2 ) denotes the parton distribution within the deutron, and in the impulse 
approximation F^ = 2F%. a(hard) denotes the cross section for Higgs meson production 
in the parton - parton collision. 

Since the deutron is a two nucleon bound state, we may also consider a production 
mechanism in which two nucleon-nucleon interactions are responsible for the Higgs pro- 
duction. These are shown in Fig. lb which, after squaring, gives the Mueller diagrams of 
Fig.lc. However, to obtain the full answer a more careful treatment is required in which 
we sum over all contributions to the inclusive production, according to the AGK cutting 
rules ||. This procedure leads to a cancelation of the contributions shown in Fig.lc, so 
that the net inclusive cross section is given by Fig. la. 

The above consideration fails to include the interference diagram shown in Fig. Id. We 
claim that this contribution survives the cancelations implied by the AGK cutting rules, 
and this results in a violation of the factorization theorem. In this note we calculate this 
contribution and assess the consequences in detail. We wish to remind the reader that the 
application of the AGK cutting rules are different when calculating inelastic and inclusive 
processes. This is shown in Fig. 2 for nucleon-nucleon scattering. 

The present paper continues the approach utilized in many cascade models for inclusive 
production in nucleus- nucleus interactions (see, for example, ||). Recently, this mecha- 
nism has been revived by Ryskin ||, who pointed out its importance for nucleus-nucleus 
collisions where its contribution is proportional to (Ai^)^, wheras using the factoriza- 
tion theorem one obtains only an (A1A2) dependence. The result obtained by Ryskin is 
evident from Fig. Id as the contribution of this diagram is proportional to (A 2 A 2 ,) - the 
number of two nucleon-nucleon collisions -divided by ttR 2 Ai ttR\ 2 - which comes from the 
integration over g x and q 2 (the momentum transfers in Fig. Id). Our paper differs from 
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Refs. |4[] ||, in as much as we have included explicitly the relevant contributions and 
cancelations implied by the AGK cutting rules. For a deutron-deutron collision, we find 
that the interference diagram resulting from the two nucleon interaction (Fig. Id) has to 
be added to the one nucleon interaction of Fig. la. Checking the numerics, we obtain a 
significant addition to the cross section for central mass production, which is in the few 
GeV range. This is relevant for estimates of the minijet cross section. 

In Fig. 3 we show the diagrams illustrating the new mechanism for a centrally produced 
Higgs meson. Fig.3a shows the process where the Higgs, which is produced centrally, is 
seperated by two large rapidity gaps (LRG) from the small final state multiplicities which 
occur on the edges of the rapidity plot. We define this process as a double Pomeron 
exchange reaction, and denote its contribution as The contribution of Fig.3b is 

denoted by er 2 ■ This is a mixed diagram, where the Higgs which is produced in a double 
Pomeron process, is superimposed on a normal uniform rapidity distribution typical of an 
inelastic nucleon-nucleon reaction. Finally, Fig. 3c describes Higgs production as part of 
the nucleon-nucleon background to the rapidity distribution. We denote this contribution 
by <j 2 • To obtain our final result we need to sum over all three of the above contributions, 
noting that these are not necessarily positive. 

The nucleon-nucleon (NN) amplitude for Higgs production via double Pomeron ex- 
change, has been calculated 0, to be 

A H = A(NN -> iV + (LRG) + H+ (LRG) + N) = 2g H A P , (2) 

gn is the the vertex of the hard parton - parton — > Higgs process and Ap is the Pomeron 
exchange amplitude. Fig. 3d illustrates graphically the source of the factor 2 (in Eq. (0)), 
which plays an important role in our calculation. The second ingredient in our calculation 
is the amplitude shown in Fig.3e. This amplitude has no analog in the case of a single 
nucleon-nucleon interaction, and it depicts the cut in the diagram shown in Fig. 3b. We 
note that this diagram is equal to ImAn, unlike the case for inelastic nucleon-nucleon 
cross section where <jj n = 2ImAp. The above result follows from the unitarity constraint 

Recalling that the integration over the longitudinal components of the vector in 
Fig. 3b results in a negative sign for the interference diagram (see Ref. || for details), we 
can easily calculate all diagrams contributing to the inclusive Higgs meson production. 
Indeed, calculating af\ the sum of the diagrams in Fig. 3a, and using Eq. (^]), one has 

4 0) = 8<?| { (ReAp) 2 + (ImApf } . (3) 
For the sum of diagrams in Fig. 2b we have 

4 1] = -16g 2 H (ImA P ) 2 . (4) 
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(2) 

Finally, for , shown in Fig. 2c, we obtain 

af ] = 8g 2 H (ImA p ) 2 . (5) 

The relation between the different contributions is 

a® = (l + p 2 ):-2:l, (6) 

where p = f^ p p ■ This should be compared with the AGK relation for the inelastic 
cross sections 

*£ } = ■■ $ = (1 + P 2 ) : "4 : 2. (7) 

Summing all contributions we obtain an additional term in the cross section compare to 
the one given in Eq.(l) 

a(D + D -> H + X) = 8g 2 H (ReA P ) 2 . (8) 

We would like to draw attention to the novel fact, that the contribution of the new 
mechanism is proportional to the real part of the amplitude. For a soft Pomeron with 
ap(0) = 1, the contribution of Eq.(8) vanishes, and we recover the factorization theorem 
[IJ, i.e. an exact cancellation of the diagrams shown in Fig.lc. 

The detailed kinematics of a double Pomeron Higgs meson production, as well as its 
amplitude, have been discussed in Ref.|J. The most important kinematical feature is 
that at high energy and for a Higgs meson in the central rapidity region: y^> 1, and Y - 
y 3> 1 ( see Fig.4a). Accordingly, the momentum transfers q\ and q2 are a small fraction 
of the longitudinal components, j|4 and |M 1. We can, therefore, write q\ = q\ t 

and q\ = q\ t . Expanding these momenta with respect to pi and p 2 we have 

qin = xiPin + /3iP2/i + qu;n , (9) 

qiH = X 2 P2fi + «2Pl At + <?2t; M • 

We obtain 



<M* + ( qi -q 2 y t 
xi = 1 F e y > a 2 , (10) 



M^ + {q L -q 2 ) 2 t 



x 2 = 1 -= e v > fa , 



where M denotes the Higgs mass. 
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We write the expression for the diagram of Fig. 4a assuming that k t ^> qu ( see 
Fig. 4a). Indeed, q u ~ q2t ~ R , where R is the hadron (nucleon or nucleus) radius, while 
k t can be as large as M. Inserting the vertex for Higgs emission by a gluon - gluon fusion 
0, namely T H = gH(ki^k 2u — g^ v kx-k 2 ), we obtain the following formula for the Higgs 
meson production amplitude of Fig. 4a 

Ah = \g H J dk%(j)(xi, kl) 0(x 2 , k 2 ) F x {q\AF 2 {q\A , (11) 

where a s (k 2 )xG(x } k 2 ) = J k dk l2 a s (k' 2 )<p(x, k' 2 ), F 1 and F 2 are the form factors describ- 
ing the t dependence of the Pomeron - target vertex (-Fi(O) = -F 2 (0) = 1) . We wish to 
stress that (f>(x, k 2 ) are not the usual gluon densities, but the asymmetric (off diagonal) 
density functions that have been introduced in Ref. 0, and have been discussed in more 
detail in Ref.|J. In the leading log(l/x) approximation of perturbative QCD ||, they 
can be approximated in the region of very small X\ and x 2 by the usual nucleon gluon 
densities. 

Neglecting the k 2 dependence of as(k 2 ), = dxG ^ k ' ^ (& 2 ) <7>_1 , where < 7 >, the 
average anomalous dimension, can be considered to be a smooth function of both x and 
k 2 in the region of small x. This is a semi-classical approximation, which holds at small 
x H, and we can use it to estimate the integral in question. Substituting in Eq. (|TT|) 
we see that for all < 7 > < | the integral is infrared divergent. Only large distances 
contribute to the integral and we can, therefore, replace this integral by the exchange of 
a soft Pomeron as was done in Ref.||. In most of the models for the soft Pomeron [|i~0|j , 
the real part of the Pomeron amplitude is very small, therefore, we expect that the new 
mechanism will only have a very small contribution to soft inclusive production. 

In general, < 7 > is a function of x and k 2 , approaching the value < 7 >= \ in some 
kinematic region. In the vicinity of < 7 >= | the integral over k 2 becomes logarithmic 
and all distances (rj_), down to very small r± « -j, contribute to the integral. In this 



case we need to consider two new physical phenomena: the BFKL Pomeron |IT] and 
shadowing corrections (SC) [|l!|. For the BFKL Pomeron, the value of the anomalous 
dimension cannot be larger than 7 = | so we evaluate our integral considering 7 — > |. In 
this case 

7T dlmAp 7T 

ReA P = ,1 n A = ^dmAp = 2a s N c \n2 , 12 



where ap(0) = l + u) P and we have used the explicit calculation of u P = uj l JTTJ]. One can 
see that for ap(0) = 1.25, ReAp w ImAp. Therefore, we expect a considerable violation 
of the factorization theorem for the BFKL Pomeron, as well as for a DGLAP Pomeron 



with an intercept which is well above unity [ 13 . 
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An estimate of SC in the case of Higgs production via double soft Pomeron exchange 
in nucleon-nucleon scattering has been discussed in [0. For the DGLAP evolution SC 
become important in the vicinity of 7 = ~, where ImAp oc (^) UJcr || with uj cr = 2N ^ as . 
Once again, for «p(0)= 1.25, we have ReAp pa ImAp. 

We can, thus, expect a considerable violation of the factorization theorem for hard 
inclusive production, in the kinematic region where the anomalous dimension approaches 
\ and where u> p is well above 0. 

To give a numerical estimate for the additional cross section we write 

*>■" _ t*to>&> M > = a(hard)l r dk fA^fl Hx ^ )? , (13) 



dy J (16vr) 2 1 R\R\ v " Jq> * dlml/a^) 

where a(hard) denotes the hard cross section, which is the same as in the factorization 

theorem (see Eq. ([!])), Ri and R2 are the radii of the colliding hadrons (nuclei) normalized 

2 

so that: Fi(t) = 1 + jRft. For the nucleus - nucleus collision 4>i = A<Aw and R 2 = r^A'f 

4 4 

and we get a NF oc A{ A| . 



In Fig. 5 we plot the ratio of the nonfactorized contribution of Eq. (|13|) to the cross 
section from the factorization theorem of Eq. ([!]) , namely 



R\Rl Xl G x {x u ^)x 2 G,{x^) ' 1 ] 

where XiGi(xi, M 2 ) is the gluon density of the z-th hadron. In Fig. 5 we took y = 0. Notice, 
that this ratio does not depend on the hard cross section, and as such is also applicable 
to any central production process, in particular, hard minijet and jet production. In 
evaluating Eq. fll4|) we use the GRV parametrization |L3| for the gluon density. We assume 



that R\ = R\ = 5GeV~ 2 (see Ref.|Pf ). For the initial virtuality we take Ql = lGeV 2 , 
as the GRV parametrization is in agreement with the HERA data on F 2 (x,Q 2 ) for 
all Q 2 > 1 GeV 2 . 

One can see, that R is quite big for low masses and decreases when M increases. For 
example, R = 1.2 for M — 10 GeV. We do not expect the Higgs meson mass to be that 
small, but our result is also applicable in the case of jet production (see Fig. 4b) where 
M ~ 2pj_. Accordingly, with such a value of R we can expect minijet production (jets 
with p± m 5GeV ) which can be responsible for the structure of the minimum bias events 
at the Tevatron energy. 

The value of R is bigger for nucleus - nucleus interaction. To estimate this value we 

11 2 

need to multiply Eq. (|14D by a factor of A^f ^4| e //- Using the simple relation R\ = r^A^ , 
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one has Al ff = -$■( A3 — 1 ) + 1. Therefore, for a gold - gold interaction we expect that 
R is enhanced by factor 4-9. 

In calculating R we have not included a Sudakov form factor, which is believed to 
be necessary in a diagram of the type of Fig.5a, i.e. there is no emission of gluons with 



rapidities between yi and y 2 in Fig. 4 (see Ref. [16] for details). However, as we have 
calculated an inclusive production process, this type of emission should be included (see 
below). As a first approximation we calculate the ratio R without the Sudakov form 
factor and neglect the BFKL enhancement, which can be important for large values of M 
(see Fig. 6). 

The diagrams of Fig.4 are only a low order approximation for the new process. Many 
additional gluons can be emitted, which lead to more complicated diagrams, shown in 
Fig. 6. The emission of the additional gluons can be summed, and the cross section 
can be described by a Mueller diagram with the Pomeron - Pomeron interaction (A) as 
seen in Fig. 6. It is important to note that we get the same equation (see Eq. (|8|)) for 
the inclusive Higgs cross section, when using the AGK cutting rules for the Pomeron - 



Pomeron amplitude (A) as was obtained in Ref.JlTj. Indeed, it was shown in Ref.fT7f that 
for different cuts of A we have a transition matrix which can be written in the following 
form 

XI: X\ : XI: XI = I: I: 1 -: I (15) 

with all other X{ are equal to zero. Here, we denote by X{ the elements of the transition 
matrix for the Pomeron - Pomeron amplitude, where i (f) is the initial (final) multiplicity 
state (say, above (below) the Pomeron - Pomeron interaction in Fig. 6) and i (f) can be 
equal to 0,1 and 2, respectively. Using this transition matrix together with the AGK 
cutting rules for the inelastic cross section, one obtains the inclusive cross section for 
Higgs production 

-j-{D + D — ► H + X) (16) 
ay 

= p 2 ^j-{D + D -> Mx{n « 1) + (LRG) + (H + X) + (LRG) + M 2 (n « 1)) , 
dy 

where Mi and M 2 denote the bunches of final particles with low multiplicities. Eq. (|i~6|) 
provides a method for calculating the inclusive Higgs production in terms of the cross 
section for the process with two large rapidity gaps, and the Higgs meson in the central 
rapidity region accompanied by a bunch of final hadrons. 

The main result of our paper is Eq. (|T4]), which gives the scale for the violation of the 
factorization theorem due to the new mechanism present in the inclusive production. It is 
clear, that this formula describes not only the process of a heavy Higgs meson production, 
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but also other central hard processes. The only difference is in the explicit expression for 
the hard cross section in Eq. fll3|) . From Eq. (|l~4Tj one can see that the new mechanism 
only gives a small contribution to the inclusive cross section for the production of light 
hadrons with small transverse momemtum (soft inclusive production), because the real 
part of the elastic amplitude is very small at high energy. For hard inclusive production, 
the new mechanism yields a sizeable contribution, especially in the kinematic region where 
the anomalous dimension of the gluon density is close to |. 

Our general conclusion is that the factorization theorem]]]] is only approximate. We 
suggest that its proof should be reexamined with the view of finding the source of the 
violation. 
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Figure Captions 



Fig.l: Mueller diagrams for inclusive Higgs production in a deutron - deutron collision : 
(a) the factorization theorem contribution, (b) the two nucleon interactions, 
(c) the diagrams, that cancel in total inclusive cross section, (d) the interference 
diagram that should be addded to factorization theorem contribution. 

Fig. 2: The AGK cutting rules for inelastic and inclusive nucleon - nucleon cross sections. 

Fig. 3: Mueller diagrams for the new mechanism of the inclusive production: 
(a) , (b) (c) and the relation between different 
Mueller amplitudes ((d) and (e) ). 

Fig. 4: Feyman diagrams (a) for Higgs and (b) for two high px jet production. 

Fig. 5: Calculated R (see Eq.(14)). 

Fig. 6: The general Mueller diagram for the new mechanism 
for the hard inclusive production. 
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Fig. 7 



